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Abstract 

We calculate the 0{a^g) massive operator matrix elements for the twist-2 operators, which 
contribute to the heavy flavor Wilson coefficients in unpolarized deeply inelastic scattering 
in the region ^ m?. The calculation has been performed using light-cone expansion 
techniques. We confirm an earlier result obtained in [1]. The calculation is carried out 
without using the integration-by-parts method and in Mellin space using harmonic sums, 
which lead to a significant compactification of the analytic results derived previously. The 
results allow to determine the heavy flavor Wilson coefficients for F2{x,Q'^) to O(a^) and 
for Fl{x,Q'^) to O(a^) for all but the power suppressed terms oc {m? /Q"^)^ , k >1. 



*Dedicated to the Memory of W.L. van Neerven 



1 Introduction 



Deeply inelastic electron-nucleon scattering at large momentum transfer allows to measure the 
parton distribution functions of the nucleons together with the QCD scale Aqcd- In the region 
of large hadronic masses ~ Q"^!^! — x)/x the sea-quark distribution receives substantial 
contributions due to heavy flavor (charm and beauty) pair production. At the level of leading 
twist T — 2 their contribution to the deeply inelastic structure functions is described by heavy 
quark Wilson coefficients, which are convoluted with the light quark and gluon parton densities. 
Depending on the range of the Bjorken variable x and the gauge boson virtuality Q"^, these 
contributions can amount to 20-40% of the structure functions [2]. The unpolarized heavy flavor 
Wilson coefficients were calculated at leading order (LO) in Refs. [3]. The next-to-leading order 
(NLO) corrections were derived in semi-analytic form in Refs. [4] in x-space for the structure 
functions ^2(2;, Q^) and Q^). A fast numerical implementation in Mellin space was given 

in [5]. For the asymptotic region ^ m^, an analytic result for the heavy flavor coefficient 
functions was calculated to 0{al) in [1]. In the case of the structure function Fl{x,Q^), the 
asymptotic result to 0(q;^) was derived in [6]. The leading order small-x resummation for 

F2l{x,Q'^) was calculated in [7]. The heavy quark Wilson coefficients differ significantly from 
those of the light quarks even in the asymptotic region. Therefore, the scahng violations of the 
heavy flavor part in F2^l{x,Q'^) are different from those of the light flavor contributions. Both 
for the measurement of the QCD scale Aqcd and for the extraction of the light parton densities a 
correct description of the heavy flavor contributions is therefore required. As shown in Ref . [1] , in 
case of the structure function ^2(2;, Q^) the asymptotic heavy flavor terms describe the complete 
contributions very well already for scales Q^^SOGeV^, whereas for Fl{x,Q'^) this apphes only 
at much higher scales, (5^>800GeV^. 

In the present paper, we recalculate for the flrst time the asymptotic 2-loop corrections to the 
heavy flavor structure functions using a different method than in Ref. [1] . All logarithmic terms 
and the constant term of the heavy flavor Wilson coefficients are obtained due to a factorization of 
this quantity into the massive operator matrix elements and the light parton Wilson coefficients, 
which are known from the literature [8-10]. In [1] the massive operator matrix elements were 
derived in momentum-fraction (a;)-space. The corresponding 2-loop integrals were simplified 
using the integration-by-parts method [11]. We will work in Mellin-space, accounting for the 
appropriate symmetry of the problem, and do thoroughly avoid the integration-by-parts method. 
This requires to solve more complicated Feynman-parameter integrals. However, in this way we 
are able to derive by far more compact results, even for the individual Feynman diagrams. In the 
direct calculation, we use Mellin-Barnes integrals [12,13] and representations through generalized 
hypergeometric functions [14]. A brief account on scalar 2-loop integrals to be derived in the 
present calculation was given in [15-17]. The final expressions obtained allow to represent the 
asymptotic heavy flavor contributions to the deep-inelastic structure functions in Mellin space in 
a completely analytic form. Precise representations of the analytic continuations of the harmonic 
sums w.r.t. the Mellin index N to complex variables are given in [18,19]. A fast single numerical 
inverse Mellin transformation yields the structure functions in x-space. This representation is 
well suited for fast numerical data analysis [20]. 

The paper is organized as follows. In Section 2, we describe the principal method applied 
to derive the 2-loop corrections in the asymptotic region ^ rn?, covering all contributions 
but the power corrections oc {rn^ /Q'^Y, k > 1. In Section 3, the renormalization of the massive 
operator matrix elements is described. Section 4 gives a brief account of the 1-loop corrections. 
The 2-loop corrections to the operator matrix elements are derived in Section 5. Working in 



2 



D — 4+£ space-time dimensions, the sphtting functions, related to the problem, can be unfolded 
in leading and next-to-leading order, which provides a check for the calculation. Here we also 
discuss the mathematical structure of the results and compare to the result obtained in Ref . [1] . 
Section 6 contains the conclusions. The appendices summarize details of the calculation and 
different types of summation formulae used in the present calculation, which are of general 
interest for other higher order calculations. 



2 The Method 



In the twist-2 approximation, the deep-inelastic nucleon structure functions Q^), i — 2,L, 
are described as Mellin convolutions between the parton densities fj{x,ii^) and the Wilson 
coefficients Cj (x, Q'^/ jj?') 



to all orders in perturbation theory due to the factorization theorem. Here /x^ denotes the 
factorization scale and the Mellin convolution is given by the integral 

[A®B]{x)= f dxi f dx2 6{x - X1X2) A{xi)B{x2) . (2) 
Jo Jo 

Since the distributions fj refer to massless partons, the heavy flavor effects are contained in 
the Wilson coefficients only. We will derive the massive contributions in the region ^ 
m^. These are the non-power corrections in m'^/Q'^, i.e. all logarithmic contributions and the 
constant term. We apply the coUinear parton model, i.e. the parton 4-momentum is p = xP, 
with P the nucleon momentum. The massive Wilson coefficients itself can be viewed as a 
quasi cross section in pV* scattering, where V* denotes the exchanged virtual vector boson. 
In the limit ^ m^, the massive Wilson coefficients H^'^^{Q'^/m'^,m^/iJ,'^,x) factorize [1] 
into Wilson coefficients C^i-k (Q^//^^)^) accounting for light flavors only and massive operator 
matrix elements ^^f^ {m^ / ^^,x). 

The latter take a similar role as the parton densities in (1). They are process independent but 
perturbatively calculable. The factorization (3) is a consequence of the renormalization group 
equation. The operator matrix elements ^^'f ^ obey the expansion 



(2 \ °° 
^ / 1=1 



of the twist-2 quark singlet and non-singlet operators O^^^ between partonic states which 
are related by coUinear factorization to the initial-state nucleon states \N). Gg = Q;s(/x^)/(47r) 
denotes the strong coupling constant. The Feynman rules for the operator insertions are given 
in Figure 1. Since the operator matrix elements are process-independent quantities, the process 

S NS 

dependence of -f^2,'L;i described by the associated light parton coefficient functions 

) = E <ci\ , NS, S, , . (5) 
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The MS coefficient (and splitting) functions, in the massless hmit, corresponding to the heavy 
quarks only, are denoted by 



C 



C2,L;k{^:NL + NH 



(6) 



where Nh, Nl are the number of heavy and light flavors, respectively. In the following we will 
consider the case of a single heavy quark, i.e. Nh = 1. The formalism is easily generalized to 
more than one heavy quark species. 

The massive operator matrix elements to 0(0^) allow to calculate the heavy quark Wilson 
coefficients in the asymptotic region for F2{x,Q'^) to 0{a^) [1] and for Fl{x,Q^) to O(a^) [6]. 
The general structure of the Wilson coefficients is 



'2,9 



trPS 



trNS 



ttVS 



rrNS 



9 ' 9 



Q2 rn?_ 
m2' ^2 

m2' ij? 

g2 



9 ' 9 



Q2 ^ 

9 ' 9 



2,9 \ ^2 



+ a: 



(2) ( Q 



^PS,(2) / 



A 



(1) 



PS,(2) / 



(1) / Q 



2,q I 2 



+ A 



(2) ( 

Qg I ^2 



NS 



(2) (Q 



^Q, \ ^2 



.NS,{2) ( 



A 



(1) 



S)-s(f)-s(f 



+ a. 



,(2) / ^^ 



Qg \ ^2 



^^Qg 



^(2) / 



, ^(3) 



2/^PS,{2) / 

2/=?NS,(2) / 

^-2 



A: 



PS,(2) / 1^ 

Qq 



.NS,{2) / }J_ 

qq,Q 1^2 



, ^(1) (Q^\ , ^NS,(3) / 



(7) 
(8) 
(9) 



(10) 

(11) 

.fl2) 



3 Renormalization 

The calculation is performed in D = 4 + £ dimensions in the MS-scheme. The massive operator 
matrix elements Ajj(m^//i^), with being the renormalization scale, exhibit ultraviolet diver- 
gences which have to be removed by operator-, mass- and coupling constant renormalization. 
Furthermore, it contains coUinear singularities, in the present case beginning with 2-loop order. 
We follow the notation of Ref. [1] and briefly summarize the renormalization procedure. 

The external hues of the diagrams are treated on-shell after factorization. The scale for the 
process is set by the heavy quark mass m. The Zo j^-matrix performs the renormalization of 
the composite operator turning Aij{m^ / fj,^, Us, s) into Aij{m^/ fi^, Ug, e), 



Aij{m'^/l/,as,e) = Zo,jk{as,e) ^ Akj{m^ / i/,as,e) , 



(13) 
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with 



k=0 



7^^ 



'0,kj 



(14) 
(15) 



Here, 7/^ denote the singlet anomalous dimensions, which are related to the splitting functions 



by 



Jo 

The coUinear singularities are removed by the matrices rfcj(as,£:). 



with 



(16) 



(17) 



(18) 



k=0 



To 2-loop order, one has 

r 



Sij + asSe-Pij 



(0) 



lp(0) 

2 



(0) 



kj 



^3 



2£ 



(19) 



The factorization and renormalization scales are choosen to be equal, iir = fip = IJ'- The 
spherical factor collects all terms to be removed in the MS scheme 



5, = exp{|[7s-ln(47r)]} . 



(20) 



7s is the Euler-Mascheroni constant. 

Finally, the mass and coupling constant renormalization has to be carried out. The bare 
coupling Os and bare mass m are related to the renormalized quantities by 



cis = Z^as = as{ii^) [1 + as{iJ?) ■ 5as\ + O(a^) , 
ZrrJn — agSm + 0(al) . 



m 



We choose the on-mass-shell scheme for mass renormalization. Here, 



H 



_ Z[ + Zj 
{zi + Z^fl- 



= 1 + agSe 



e e 



E 

NH=i 



m 



Nh 



e/2- 





f 


[5-4] 






e 



(21) 
(22) 



(23) 
(24) 
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in Feynman gauge [22-24]. Nh denotes the heavy quark species and the SU{3)c color factors 
are Cp = {N^ - l)/(2A^c), Ca = Tr = 1/2, = 3. The Z-factors in (23, 24) read : 



2 4 

3 3^^ 
5 4 



as — -TR 
e 3 



1 + 



E 

Nh=4: 



m 



The lowest expansion coefficients of the /3-functions are 



Nh 



e/2 



(25) 
(26) 

(27) 



/3o = -[llCA-4TRAr^], 
4 

Po,Q — ~'^'^R ■ 
The bare and renormalized coupling and quark mass are related by 



Sag = Se 



e ^ e 



m 



Nh 



e/2 



1 + 



5m — CpSrin ( 

11^ 



m 



2\ £/2 



5-4 



The renormalized operator matrix element Aij is now given by 



A. 



Sij + a 



ik 



-1,(2) 



{ 



+ 0(a3). 



(28) 
(29) 



(30) 
(31) 



(32) 



In the following sections we first calculate the un-renormalized operator matrix elements Aij 
from which A^j is derived through (32). Due to (16, 19) the following leading and next-to-leading 
order splitting functions arc needed. We will mainly work in Mellin space and therefore list these 
functions in this representation. The leading order splitting functions read [25] 



ACf 



-25i(7V-l) + 



(AT - l)(3Ar + 2) 
Z/V(7VTT) 



^TrNf 



N^ + N + 2 
N{N + 1){N + 2) 



8Ca 



-S,{N-1)- 

N^ + N + 2 
{N -1)N{N + 1) 



N^-3N-A 



{N -1)N{N + 1){N + 2) 



2Po 



(33) 
(34) 
(35) 
(36) 
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Furthermore, the following next-to- leading order splitting functions contribute [21,26] 

5N^ + 32N* + 49N^ + 38N^ + 28N + 8 



pNS,(l)/^x - pNS,(l) _ 



qq,Q 



{N -1)N^{N + iy{N + 2f 



(37) 



CfTr i ^S,{N - 1) - ^S2{N - 1) 

4 (iV - 1)(3A^ + 2)(iV2 - llA^ - 6) 

TV^ + TV + 2 4 



(38) 



A^(A^ + l)(A^ + 2) 



where 



+ 



+ 15N^ + 36A^^ + 51A^=^ + 25N'^ + 8A^ + 4 



+IQCaTr 
+4 



Ar3(7V + l)3(7V + 2) 
/V^ -I- A/" + 2 

[5i2(iV) + 52(iV)-C2-2/3'(iV + l)] 



Ar(Ar + l)(Ar + 2) 
2A^ + 3 



(iV+l)2(iV + 2)2 



^i(iV) + 



Pi(iV) 



{N -i)m{N + i)^{N + 2y 



Pi{N) = iV^ + 6Ar8 + 157V^ + 257V^ + 36Ar5 + 85Ar4 + i28iV3 + 104iV2 
+64A^ + 16 . 

Here, the harmonic sums [27, 28] are given by 



AT 



= E7 = V'(A^ + l) + 

SkiN) 



1=1 

N 



j-^l'^ {k-l)\ 



lE 



^(fe-i)(iV+l) + C,, k>2 



S-i{N) 
S-k{N) 



(-l)^/?(iV + l)-ln(2) 



(^-1)! 

Cjfc denotes the Riemann ^-f unction. 



p'^-\N+i)-(i-^^]a, k>2. 



(39) 



(40) 



(41) 

(42) 

(43) 
(44) 
(45) 



4 The one-loop massive operator matrix elements 

At one-loop order, only gluonic terms contribute to heavy flavor production. The complete 
calculation of the differential scattering cross section (f'a{'y* + N —>■ QQ) / dxdQ^ was performed 
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in Refs. [3]. In the limit » m^, the heavy flavor Wilson coefficients are in the MS-scheme 



H 



(1) 



Q 



C 



(1) 



H. 



(1) / Q 



2,9 



9 ' 9 



m 



+ a 



2,9 



(46) 
(47) 



Here C(2 l) denote the massless one- loop Wilson coefficients and y4Q^(m^//x^) is the 

massive one-loop operator matrix element. Since Fl{x,Q^) is coUinear finite at leading order, 
the Wilson coefficient (^^^(Q^/At^) is universal and independent of the scheme or method of 

calculation. if|^^((5^/m^, rri^/ /i^) does therefore not contain contributions due to operator matrix 
elements at this order. The logarithmic contributions to Fl{x, Q'^) emerge as (m^/Q^)-ln(Q^/m^) 
and vanish in the limit ^ m^. The massive operator matrix element ^^^(m^/ Q^) is calculated 
from the diagrams in Figure 2. The symbol ® denotes the operator insertion, cf. Figure 1. The 
massive operator matrix elements are obtained by contracting the diagrams with the projector 

,2 



m 



Aqs I s, — 



1 



1 



N^-lD-2 



(48) 



where a and h are the outer color indices and /i and u are the Lorentz-indices and Gij is the 
respective Green's function. The diagrams yield 



l + (-l) 



N 



e/2 



m 



1 



{2 + e)e 




X 



2{N^ + 3A^ + 2) + £{N^ + N + 2) 



N{N + 1){N + 2) 



l + (-l)^ 



e/2 



m 



exp 



(2 + s)s " " [ 




(A^ + l)(A^ + 2) 



(49) 
(50) 



For the unpolarized operator matrix elements only the even moments contribute due to the 
current crossing relations, see e.g. [29]. The analytic continuation is performed starting from the 
even moments. The matrix elements have to be expanded to 0{e), since these terms are needed 
for the 2-loop result: 



^Q9 



1 

a. 



— —SsTr 




eVl S{N^ + N + 2) 
2) j N{N + l){N + 2) 



{N^ + N + 2) 
N{N +1){N + 2) 



(51) 



Note that the term 0{e^) vanishes. In x-space the result reads 



e/2 



"■Qg 



1^^ 



(52) 
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where 



''Qg 



(53) 
(54) 



and Pqg^ (,t) denotes the leading order sphtting function for the transition g ^ q for one (heavy) 
flavor. At one-loop order, the renormalized operator matrix element ^4^^ is obtained by 

id) ^ 7-1.(1) ^55^ 



m 



where 



7(1) 

'0,Qg 



.1p(o) 



(56) 



Here, removes the ultraviolet singularities. At this order no collinear singularities are 
present. Choosing /x^ — the massless Wilson 1-loop coefficients in the MS-scheme [8] are 
given by : 



In 



\-z 



ct:g{^) = i6r«^(i-^). 

The asymptotic heavy quark Wilson coefficients (46,47) are 



Q 



m 



Q 



(57) 
(58) 

(59) 



In 



I- z 



+ ln ( ^ 



- 32^2 + 32^-4 



(60) 



These are precisely the expressions obtained in the hmit rr? jQ"^ — > from [3] up to the constant 
term {m^jQ'^f : 



E 



compl. 



,(1)L ^ 



^ = 16T, 



R 



m 



vz(l - z) +2—2;^ In 



1 -V 

1 + v 



(61) 



h: 



compl, 
2,ff 



1 . .2 -^"^^ 



__+4z_4z^_2— ^(1-z) 



+ 



In 



1 -t; 
1 + v 



, (62) 



where denotes the cms- velocity of the heavy quarks. 



v^i 1 



Am? z 
l-z 



(63) 
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5 The two-loop massive operator matrix elements 



There arc three classes of 2-loop contributions to the massive operator matrix elements: the 
gluonic contributions (diagrams Figure 3), the pure-singlet contributions (diagrams Figure 4) and 
the non-singlet contributions (diagrams Figure 5). The diagrams are either one-loop insertions 
into one-loop diagrams or, in case of the gluonic contributions, also genuine two-loop diagrams. 
The calculation is performed using FORM [30] and Maple procedures. We calculate all diagrams 
directly, i.e. without decomposing them using the integration-by-parts method as done in [1]. 
The integrals to be performed are more involved. However, we avoid a large proliferation of 
terms in the results, in our case of nested sums of different kind, which add up to zero. Even 
in the case of individual diagrams, which are calculated in Feynman gauge, only a very small 
number of harmonic sums contributes finally. One of these, S2,i{N), does not emerge in the 
operator matrix elements. The results for the individual non-renormalized diagrams are given in 
Appendix A. It turns out. that some of the diagrams can be easily calculated to all orders in e. 
The diagrams can be represented in terms of linear combinations of generalized hypergeometric 
functions [14]. Due to the given topologies the most complex function is 3^2(01, a2, 03; bi, 62; 1)- 
In these representations, the conformal mapping of Feynman parameters is essential. The scalar 
integrals associated to the genuine two-loop diagrams have also been calculated using the Mellin- 
Barnes technique [12, 13], cf. Refs. [15, 16], and checked for fixed moments using the package 
MB [31]. The ^-expansion can be performed prior to the summation. It results into finite and 
infinite sums of various types, including harmonic sums attached with Euler Beta-functions and 
binomials. Here, we face a more general situation than in massless calculations, despite the fact 
that we work in the limit ^ m^. The sums are given in Appendix B. They were performed 
using suitable integral representations or using difference equations. Finally, the results depend 
only on harmonic sums, which are reduced further applying their algebraic relations [32]. 

After mass renormalization, the operator matrix element /iQ^^m^//x^,£^ is given by 



''Qg 



N„=4 



m 



2\ e 



N^ + N + 2 



R^F- 



iV(iV + 1)(A^ + 2) 
+ iv + 2 



-ASi{N) + 



3Ar2 + 3jv + 2 
N{N + 1) 



+ -{4:TrCf 



A^(7V+ l)(Ar + 2) 
iV2 + iV + 2 



N'^ + N + 1 



{N -1)N{N +1){N + 2) 



S2{N)-Sl{N)'^ +4- 



A^(A^ + l)(A^ + 2) 
+ 15N^ + 36A^^ + 51N^ + 25N'^ + SN + A 



7V2 



A^3(7V + l)3(7V + 2) 



+TrC, 



+ N + 2 

) I I 

' N{N + 1){N + 2) 



-2(5\N + 1) + S2{N) + Sl{N) - C2) 
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-32 



2A^ + 3 



(A^+ l)2(iV + 2)2 ' (A^ - l)iV3(iV + 1)3(2 + iV)3 



5i(Ar) -8- 

C2 



A(iv) 



£/2 



(1) 
Q9 • 



(64) 



For the ghost contributions in Figure 3, the projector reads 



m-\D-2 



— (A.p)~^ (^a6,{2),ghost _ 



(65) 



From the terms oc 1/e^, l/e, one may determine the respective QCD sphtting functions to two- 
loop order (33-39) , which are recalculated in this way. 

The constant term in 74Q^^m^//x^,£^ reads : 



iV2 + iV + 2 



TV (TV + 1) (TV + 2) 
-S^{N - \)S^{N - 1) - 2(2^1 (iV - 1) 



+ 



Ar4 + i6A^3_^15;Y2_g^_4 
H 9 S2{N - 1 

(AT + 1)2 (AT + 2) ^ ' 

3Ar4 + 2Ar3^3jY2_4^^4 

+ -^T-TT 7^ C2 



2Ar2 (iV + 1)^ (iV + 2) 

- iV3 - l6iV2 + 2iV + 4 



+ 



(AT + 1)^ (AT + 2) 



5i(Ar-l) + 



A(iV) 



+4r T / ^' + ^ + 2 



4M 



Li2(x) 
1 + x 



2Ar4 (AT + 1)^ (AT + 2) 
1 



(AT + 1) + -Sl{N) + 352(iV)5i(A^) 



+ 3'53(iV) + /9"(A^ + 1) - 4/3'(Ar + 1)^1 (AT) - 4/3(Ar + 1)C2 + Ca 



A^=^ + 8A^2 + llA^ + 2^2 
A^(A^+ l)^(A^ + 2)^ ' 



IN^ + 21A^4 + 13A^3 + 21A^2 ^ + 16 



52(iV) 



(A^-l)Ar2 (AT + 1)2 (AT + 2)^ 
A^e + 8 A^5 + 23 A^^ + 54 A^3 ^ 94 ^2 ^ 72 AT + 8 

Q Q '-'l ( A/ I 

A^ (A^ + 1)^ (AT + 2)^ 

(Ar2-Ar-4) ^(iV) . 

(A^ + 1)' (AT + 2)'^ ^ ^ (A^ - l)iV4(Ar + l)4(Ar + 2)4^" ^''''^ 
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The polynomials in Eqs. (64, 66) read 



P^(N) = N^+6N^+15N'^+25N^+3QN^+85N'^+128N^ + 10AN^ + 6AN + 16 , (67) 

p2iN) = 12N^ + 5AN^ + 136A^^ + 218N^ + 221A^^ + llOiV^ - SN"^ - 2AN - A , (68) 
p3{N) = 2N^^ + 20 N^^ + 86 + 192 + 199 - - 297 - 495 A^^ 

-514Ar^-488 7V2-416 7V2-176 7V- 32 . (69) 



The pure-singlet A^^'^'^^ yw? / ix^^ej operator matrix element is given by 

^(^■^ - -Kf)l^{-5A>-i?}4{-5^=''"}- 

16{N^ + N + 2f 



PS, (2) 
Qq 



Si 



m 



— ^TrCf 



{N -1)N^{N + lf{N + 2) 



— -TrCf 

e 



S{5N^ + 32A^4 + A9N^ + SSN"^ + 2SN + 8) ps, 



{N-l)N%N + lf(N + 2)^ 
In this case and the non-singlet case, the projector 



+ a 



Qq 



(2) 



(AT) . (70) 



A 



(2) 

Qq;{qq,Q) 



m 



Nc 4 ^ ^ r 



(9) 



(71) 



is apphed. Here i and j denote the matrix-elements of the Gell-Mann matrices. The constant 
term is obtained by 



PS,(2) 



(AT) 



TrCf< -a 



+ 



{N^ + N + 2f 



(TV - l)iV2(Ar + l)2(iV + 2) 
4P4(iV) 



{2S2{N) + C2) 



{N - 1)N^{N + 1)^{N + 2)3 J ' 

P^(N) = A^^° + 8A^^ + 29A^^ + 49A^^-11A^^-131A^^-161A^^ 
-160A^^ - 168A^^ - 80A^ - 16 . 



(72) 
(73) 



Finally, the non-singlet operator matrix element A^^'Q^ym'^/ij,'^,sj reads 

^■(S)*[?{-''-<'}-j{-i«'}+"Sl 

3A^2 + SN + 2 



12- 



+-TrCf 
e 



32 



^^ar.r^ a r .r^ ^^N" + + A7N^ + 20N - 12 



+a 



NS,(2) 

qq,Q 



(N) 



(74) 
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The constant term is given by 

%g,Q (N) = CfTr'> --S3{N) - -C2Si{N) + —S2iN) + 2 3^^^^-^^ C2 - -^Si{N) 
219N^ + 657 + 1193Ar^ + 763A^=^ - AON"^ - ASN + 72 1 

These results obtained in MeUin-space agree with those given in Ref. [1] in x-space, cf. [6]. 

The method apphed here allowed to compactify the representation for the heavy flavor matrix 
elements and Wilson coefficients. As shown in Appendix A the individual Feynman diagrams 
depend on the harmonic sums Si{N), S2iN), SsiN), S-2{N), S-3{N), S2,i{N), 5_2,i(iV) only. In 
the final result the sum 5*2,1 (A^) drops out. The x-space representation in [1] contains the 
following 48 functions : 



oyi — X ) 


1 
i 




m yx ) 


m yx ) 


ln(l — x) 


ln^{l-x) 


ln^{l - x) 


\n{x) ln(l — x) 


hi{x) ln^(l — x) 


ln2(x)ln(l -x) 


ln(l + x) 


\n(x) ln(l + x) 


ln^(x)ln(l + x) 


Liafl -x) 


ln(a;)Li2(l — x) 


ln(l -a;)Li2(l -x) 


Us{l - x) 


Sl,2{l-X) 




1 


1 


ln(a;) 


ln\x) 


ln\x) 


1-x 


1 + x 


1-x 


1-x 


1-x 


ln(a:) 


ln\x) 


ln\x) 


ln(l + x) 


ln{x) ln(l + x) 


1 + x 


1 + x 


1 + x 


1 + x 


1 + x 


ln(:r) hi^(l + x) 


\n^{.v) hi(l + x) 


ln(,r) hi(l - x) 


ln(,r)hr(l - x) 


ln(l - ,/-)Li2(x) 


1 + x 


1 + x 


1-x 


1-x 


1-x 


Li2(l - x) 
1-x 


ln(a;)Li2(l — x) 
1-x 


ln(a;)Li2(l — x) 
1 + x 


ln(l + a;)Li2(— a;) 
1 + x 


ln(l + a;)Li2(— x) 


U2{-X) 


U2{^X) 


ln(a;)Li2(— x) 


Lk{l - x) 


^hi-x) 


1 + x 


1 + x 


1 — X 


1 + x 


^1,2(1 - X) 


Sl,2{l - X) 








1-X 


1 + X 


1 + x 







As shown in [27], various of these functions have Mellin transforms which contain triple sums, 
which do not occur in our approach even on the level of individual diagrams. 

In the MeUin-space representation, the sums listed in Table 1 contribute to the result of the 
individual diagrams. Note, that we express single harmonic sums with negative index in terms of 
/9-functions and their derivatives, cf. [27]. They can be traced back to the single non-alternating 
harmonic sums, allowing for half-integer arguments. Therefore, all single harmonic sums form 
an equivalence class represented by Si{N), from which through differentiation and half- integer 
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relations the other single harmonic sums are easily derived. Further the equality, 



M 



LisOr) 
1 + x 



(7V + l)-C2/3(7V + l) = (-l) 



N+l 



S-2AN) + IC: 



(76) 



holds. ^ Therefore, the operator matrix element Aq^^ depends on one non-trivial basic function 
only [27]. The absence of harmonic sums containing {—1} as index was noted before for all 
other classes of (space- and time-like) anomalous dimensions and Wilson coefficients, including 
those for other hard processes having been calculated so far, cf. [19,33]. This can be seen 
if one represents the respective expressions in form of weighted harmonic sums, following an 
earlier suggestion of one of the authors. Linear representations do not allow this since they are 
non-minimal and contain algebraic redundancies. 



Diagram 


^1 


^2 


^3 


S-2 


S-3 


5'2,1 


S-2,1 


# x-space fct. 


A 




+ 












8 


B 


+ 


+ 


+ 






+ 




10 


C 




+ 












4 


D 


+ 


+ 












5 


E 


+ 


+ 












9 


F 


+ 


+ 


+ 






+ 




24 


G 


+ 


+ 












6 


H 


+ 


+ 












7 


1 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


20 


J 




+ 












7 


K 




+ 












7 


L 


+ 


+ 


+ 






+ 




13 


M 




+ 












7 


N 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


38 





+ 


+ 


+ 






+ 




13 


P 


+ 


+ 


+ 






+ 




14 


s 




+ 












7 


T 




+ 












7 


PS. 




+ 














PSfe 




+ 












7 


NS„ 


















NSb 


+ 


+ 


+ 










5 


E 


+ 


+ 


+ 


+ 


+ 




+ 


48 



Table 1: Harmonic sums contributing to the individual diagrams compared to the number of functions 
in x-space, Ref. [1]. 

The expressions for the renormalized two-loop operator matrix elements (32) are given by. 



^We correct a typo in [6]. The argument of the MeUin-transform in (66) reads + 1, not N. 
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ct. [1], 

-in'»>(iV)lng)+ag(7V) 

+a'^l(N} {2;3„ + [P™(iV) - P<«(A')] } (77) 

WP'? W} In^ ) + {-i/',7*'(iV)} In (^) 

+al'f'(N)-a'il{N)P<°\N) (78) 

+«S"('V) + 3*.oC2-p4°'W ■ (79) 

Besides the splitting functions up to next-to- leading order, the constant terms af^{N) and OLQg{N) 
determine the massive operator matrix elements (77-79). 

The asymptotic heavy flavor Wilson coefficients H2^l{N,Q'^) are then given by (7-12). The 
corresponding expressions in x-space are obtained applying the inverse Mellin transform. An- 
alytic continuations of the corresponding basic functions to complex values of arc given at 
high precision in [18,19]. The inverse Mellin transform to obtain the respective contributions 
for the structure functions is performed by a single precise numeric contour integral around the 
singularities of the problem, after convoluting with the evolved parton densities [20] . 

6 Conclusions 

We calculated the unpolarized massive 2-loop operator matrix elements, which are used to 
express the heavy flavor Wilson coefficients in the asymptotic region 3> for F2{x,Q'^) to 
0{a1) and for Fl{x,Q'^) to O(a^). We confirm the results obtained in Ref. [1]. The method 
applied in the present paper is widely different from the one used in [1]. We calculated the 
Feynman diagrams without applying the integration-by-parts method and worked in Mellin- 
space, to obey the natural symmetry of the problem. The calculation refers to nested sums in 
the first place, while in [1] the Feynman-paramctcr integrals were mapped to a single Mellin 
transform successively integrating Nielsen-type integrals. Furthermore we applied the algebraic 
relations between the harmonic sums to simplify the expressions further. The representation 
obtained for the individual Feynman diagrams was much more compact. Only a few harmonic 
sums contribute, which furthermore can be grouped into only two equivalence classes. This is 
to be compared to 48 functions in a;-space, which were needed to express the result in [1]. The 
present problem exhibits a more involved nesting if compared to massless two-loop calculations, 
since the heavy quark mass connects Feynman parameters, although wc work in the limit ^ 
m^. The representation of the Feynman-parameter integrals of the loop-diagrams in terms of 
higher transcendental functions, here generalized hypergeometric functions, before carrying out 
the £-cxpansion, proved to be essential for the compactification. In the present calculation new 
types of finite and infinite sums beyond the case of multiple harmonic sums had to be performed. 
The final results could again be expressed by harmonic sums. 



.(2) / ^ 



N 



/lPS,(2) 



^ N 



qq,Q 
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7 Appendix A: Results for the Individual 2— loop Dia- 
grams 

In the following, wc list the results for the individTial Feynman diagrams, in some cases to all 
orders in e, to demonstrate the simplicity of their structure as obtained by the present method 
of direct calculation. An overall factor dlS^{'m^ / ii^Y has been taken out. Here /xq denotes the 
initial scale to define the strong coupling constant. 

The individual contributions to the operator matrix elements for the diagrams of Figure 3 are : 



Pi (TV) 



-TT 



exp 



r(iv - e/2)r(iv) 



BiN,3) 



1=2 



sin((l + £/2)7r 

{im{N + lf{N + 2f 
+8(A^ + 1)(A^ + 2)(3Ar3 _ _ 6Ar - A)e 
+m{m^ + 12N^ - 9N'^ - 2SN - 20)e^ 
+{10N^ + SN^ + 6N^ + 24:N^ + 72N + 64)£^ 
+(2A^^ - lOA^^ - 36A^^ - 24:N^ + 24:N + 16)£^ 
+(-4A^^ - AN^ + 2N^ + 2N- 12)e^ 

+{2N^ + AN^ + 2N- 4)£^) 

1 16 8 2iV3-A^-2 



r{N + 2 + e/2)T{N + 3 - e) £(£ + 2) 



TrCf' 



+ -- 



+ 



£2Ar2(Ar+l) eN^{N +lf{N + 2) N^{N + 1) 
4Pi(iV) 



:C2 



+ 0{e) , 



7N^ + 18N^ + ISN"^ - 3N^ - 21iV2 _ _ 4 _ 



(7.1) 



TrCf< 



+ 



1 



24,52(A^) - 8Sf(N) 



N 



A^(A^ + l)(A^ + 2) 



A^(A^ + l)(A^ + 2) 

AT^ + 7A^ + 2 



16 40 
-^S,,,{N) + —Ss{N) 



-—Si(N)S2(N) - —SHN) - —SAN)C2 + 4- , 

n J 2y J 3^ u ^ jy n N{N + 1){N + 2] 



^ N''-9N + 2 8^ 
+ V(Ar + l)(Ar + 2) ^-(^^ + ]V^- 



16™^i±|^5.(TV) 



N^(N + 2y 



+32 



+ lON"^ + 30N^ + 37N'^ + 18N + 4 
N{N + lf{N + 2y 



(7.2) 



^ ^ . 18 1 4(13iV^ + 82Ar3 

^^^^^ "i^iV + £ 



82iV^ + N-6) 20 _ 2_ 

m{N+l){N + 2){N + 3) N^^ ' 7V^' 



2P2(A^) 



Ar3(7V+ l)2(Ar + 2)2(7V + 3) 



(7.3) 



16 



16N^ + 176N^ + 520N^ + 600N^ + 257N^ + 7N^ + 16N + 12 



TnC 



1 16 1 

£2 AT ^ £ 



8_ AT^ + lOAT^ + 59Ar + 42 

' N^^ ' ^ N{N + 1){N + 2){N + 3) 



2 
TV 



52(iV) + -S?(iV) 



8P3(A^) 



A^(A^ + l)3(Ar + 2)2(Ar + 3) 
+ lOA^^ + 99A^^ + 350A^^ + 486A^^ + 21 AN + 60 



(7.4) 



R 



-8 



2 



1 16(Ar + 3) 1 
e 



- 2 



-2 



+4 



£2 (AT +1)2 

3Ar3 ^ ^ ;^2A^ + 4 
Ar(7V + l)3(iV + 2) 

Ar3 _ - 8A^ - 36 
N{N +1){N + 2){N + ^) 
AN^ + 19A^4 + 31A^3 _ 307V2 



8(A^ + 2) 



Si{N) 



N{N + 1) 

9N'^ + AON^ + 71N^ - UN - 36 
N{N + 1)^{N + 2){N + 3) 



S2{N) 



Sl{N)+A 



(iV + 3) 



C2 



(iV+l)2 

'''^-''-S.iN) 



N'^{N + 1)^(N + 2)(N + 3) 
4P4(iV) 



Ar2(Ar+ l)4(7V + 2)2(Ar + 3) 
16A^^ + lllA^^ + 342A^^ + 561A^^ + 536A^^ + 354A^2 + 152A^ + 24 . 



(7.5) 



R 



+ - 



16 



Ca 



64 



{N + l){N + 2) 
5N + 2 



Ar'^^(^) + %2(jv + i)(Ar + 2) 



^i(A^) - 
Si{N) 



64 



(A^ + l)(A^ + 2) 
32 



16 ^ , 8 ^ , 16 
+^^.i(iV)-^^3(iV) + ^^^^^^^^ 



(Ar + l)2(Ar + 2) 
Si{N)C2 



+4 



(9A^ + 2)(2A^-3) 
m{N + 1){N + 2) 
16 



S2{N)+A 



2N^ -3N + 2 
m{N +1){N + 2) 



+ 16 



(Ar + l)(A^ + 2)"'^ ^Ar(Ar + l)2(Ar + 2) 
2iV3 + 12A^2 ^ 23N + 18 



TrCp 



(A^ + l)3(Ar + 2)2 
1 32 1 

£2 (iV+ l)(iV + 2) ^ £ 



^ 17iV2 + 47iV + 28 
' (AT + 1)2 (AT + 2)2 



(Ar + l)(Ar + 2) 
38 



S^{N) 



(Ar + i)(Ar + 2) 



52(A^) + 



(Ar + i)(Ar + 2) 



Sl{N) 



(7.6) 



17 



C2 - 4 ,,,,, , ,,,,,, , ^i(A^) 



+36 



(N + 1)(N + 2y^ N(N + iy(N + 2y 



-8 



(7V + l)3(Ar + 2)3 J 
CaI f 1 32 



(7.7) 



£2 (iV + l)(A^ + 2) 



16 



N + 3 



N{N +1){N + 2) 



Si{N) 



N^ + 7N + 8 



-4 



- ISiV + 9 



N{N + 1){N + 2){N + 3) 



-4 



+4 



sm - 



(iV + l)2(iV + 2)2 

_ 2iV + 9 
Ar(iV+ l)(iV + 2)(iV + 3) 
3iV4 + iV3 _ 27iv2 _ 35JY _ 34 

Ar(iV + l)2(iV + 2)2(iV + 3) 
{N - 1){UN^ + UOiV^ + 337A^2 ^ 453^ + 240) 
{N + lf{N + 2y{N + 3) 

16 ^ 16(iV + 4) 



(7V+l)(iV + 2) 
^i(A^) 



52(A^) 



C2 



(7.8) 



TrCaI ^ 



(Ar + i)(Ar + 2) 



^i(iV) 



(Ar + i)(Ar + 2)2 



1 

+ - 

s 



32 



+ 



A{AN + 3) 



{N+l){N + 2) 



{N+l){N + 2) 



+ 9N^ + 17N + 8^ 2N^ + 8N'^ + 19N + 16 

+8 . ...... . 5i(A^) - 8 



N{N + iy{N + 2y 
8{2N + 1) 



16 



(Ar + l)2(7V + 2)3 

4(18iV + 17) 



Ar + 2 



5-3 (iV) 



3(Ar+i)(Ar + 2) 



{N + l){N + 2) 
+J^S_2{N)S^{N) + ^M±!) 52(A^)5i(A^) - 



(A^ + 2) 

32 
~ N + 2 

5-3 (iV) 



S-2{N) 



5-2,1 (iV) 



Ar + 2 



(iV+l)(A^ + 2) 



3(Ar + l)(Ar + 2) 



+ 



(7V+l)(7V + 2) 



C25l(iV) - 



16(iV2 -N-A) 
(iV + l)(Ar + 2)2 



5_2(iV) 



47V4 + 7V3 _ 7^2 ^ 7^ ^ g 



52(Ar) + 2 



3A^3 + 77V2 - 3Ar - 8 



5?(iV) 



A^(A^ + l)2(Ar + 2)2 ' N{N + iy(N + 2y 

4(A^ + 4) , ,4A^5 + 36A^^ + 114A^3^174A^2^;^37^^48 



rC2 



(A^ + l)(A^ + 2)2^" N{N + lY{N + 2f 

SN^ + 68A^^ + 247 ^ 449^2 ^ 4037V + 144 1 



5i(7V) 



+4 



(Ar + l)3(7V + 2)4 
16 



7 



64 



N{N + 2) 



(iV + l)(iV + 2) 
5i(iV) + ^28 



S2{N) + 



16 



(7V + l)(7V + 2) 



(7V+l)(7V + 2) 



52(iV)5i(7V) + 



(iV + l)(iV + 2) 
32 

~ 3(iV + l)(Ar + 2) 

sm 



sm 



Ss{N) 



3{N+l){N + 2) 



18 



8(3iV + 2) 



A^(A^ + l)(A^ + 2) 

16{5N^ + 9N + 6) 
' N{N + iy{N + 2) 



S2{N) - 



Si{N) - 



8{3N - 2) 
A^(A^ + l)(A^ + 2) 

192 

{N + l){N + 2) 



S!{N) 



(7.9) 



ZIrCa exp (^Z^ ^ ^ J ^ 2 + £/2)r(iV + 3-5) e{e + 2) 
(|4(7V + 2)(4Ar2 + 4Ar _ 5) _ 4(llAr2 + QAT + 9)^ - {AN^ - 2N^ - 27N - 2)e^ 
+ (4A^2 + 2A^ + 9)£^ - 2(2A^ - 1)£^) 



TrC, 



1 8(4iV2 + 4iV - 5) 1 4(4iV5 + 22Ar4 + hn^ + i3Ar2 + 35iV + 10) 



£2 A^2(7V+1)2 



+ - 

2 rT7^C2 



A^3(7V + l)3(Ar + 2) 
2P5(iV) 



4iV2 + 4iV-5 

iV2(iV + l)2 > " iV2(Ar+l)2 iV4(Ar + l)4(Ar + 2)2 

20Ar^ + 64iV^ + 120iV^ + "^AN^ - \Am^ - 2537V2 - 1007V - 20 . 



+ 0(£), (7.10) 



ATrCa exp [ ± ^.A ,^(^+^7(^)^(^-^) - 23iV - 

\^ i jT{N + 2 + e/2)T{N + 3-e) e(£ + 2) 

-{7N^ + 9N + 36)e + 2{N^ + AN + l)e'^ + {AN + 9)e^ + 2e^^ 

1 8(3A^2 _ 2SN - 20) 1 4(10Ar^ + 7A^=^ + 517V2 + 172A^ + 112) 



20) 



TrCa 



+A 



£2 (AT - l)Ar(Ar + 1)2(7V + 2) e 

3A^2_23Ar-20 ^ 

S2{N) + 2 



{N -1)N{N + lf{N + 2)^ 
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(A^- l)A^(A^ + l)2(Ar + 2) 
2P,{N) 



(N -1)N(N + 1)^(N + 2) 



C2 



+ 0{e) , 



{N - 1)N{N + iy{N + 2^ 
lAN^ + 56N^ + 153A^^ + 139A^^ - AlAN^ - 908N - 448 



(7.11) 
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16 ^ . 4iV6 + 30iV5 + 55i^ 



-Si{N)-A 



AN^ + 30N^ + 55A^^ + 38iV3 + 4iV2 - lOiV - 4 

A ra / A r I 1 \a/ AT i o\ 



A^(A^ + 1)'^^^''^ A^3(7v + l)3(Ar + 2) 
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-S2iN) - 



N^{N + 1) 

(N + 2){2N + 1) Pr{N) 

N^{N + 1)^ m{N + lY{N + 2) 
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8A^^ + 68A^^ + 16AN^ + lllN^ + 78N^ + 12N^ + UN'^ + UN + 4 . 
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N{N +iy{N + 2^ ' (Ar+ l)(iV + 2)2 

2N^ + 48A^^ + 174Ar3 + 2A2N^ + 161N + 64 

N{N + l)%N + 2f ^ 
lON^ + 92A^^ + 329A^=^ + SSIA^^ ^ 5077V + 176 1 



C2 



{N + lf{N + 2y 



+ 



1 4(2Ar3 + N^-3N-1) 



(7.16) 
(7.17) 



N^{N + 1)' 



;S2{N) 



2Pi2{N) 



N^{N + 1)^^ m{N + iy{N + 2) 

AN^ + AN^ - SN^ - 2N^ + 16A^^ + 9A^ + 2 . 

1 S{N^ + 3A^ + 4) 1 A{2N^ + 5A^^ - 3N^ - 20N - 16) 

£2 (AT - 1)N{N + iy{N + 2) ~ e (A^ - 1)A^(A^ + l)3(Ar + 2)2 



(7.18) 



TrCa< 



+4 



A^2 + 3Ar + 4 



(AT - l)Ar(Ar + l)2(Ar + 2) 



^2(A^)+2 



A^2 + 3Ar + 4 



{N -1)N{N +1)^{N + 2) 



C2 
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{N- 1)N{N + l)4(iV + 2)3 
Pi3(iV) = 2N^ + AN^ - 13N^ - 35N^ + UN^ + 92N + 6A . 

The pure-singlet contributions read : 

1 16(iV2 + iV - 2) 18(5A^3_5^2_^g^_4) 
£2 m{N+iy e N^{N +lf{N + 2) 

4Pi4(iV) 



Pl4{N) 



%{n)-a ZjJ\ , ' c(2) + 



iV2(iV + 1)^ 



■iV2(Ar + 1)2 



iV4(iV + l)4(iV + 2)^ 



iV^ + TAT^ + leAT^ - 9N^ - 26N^ + 61N^ + UON^ + UN + 8 . 



(7.19) 



(7.20) 



128 



1 



128(3 + 2N) 



£2 (N - 1)N{N + l)(iV + 2) e{N- 1)N{N + 1)2(A^ + 2)2 
64 „ 32 



-S2{N) 



{N -1)N{N +1){N + 2) ' ' {N -1)N{N +1){N + 2) 
32{N^ + 6N^ + iV2 - 2AN - 24) 



C2 



(AT - l)iV(iV+ l)3(iV + 2)3 



Finally the non-singlet contributions are : 



(7.21) 



18{N^-2 + N) 1 8{N^ + 2N^ - WN^ - 5N + 3) 



£2 3N{N + 1) 



9N^{N + iy 



IIN^ + 33N^ -34N^ -57N^ + 5N^ + 6N -9 N^-2 + N ■ 
—4 -TT 2 —Q2 ) ■ (7.22) 



27Ar3(iV+ l)s 



3iV(Ar+ 1) 



32^,,,, 32 



1 

+ - 

e 



16^,,,, 80^,,,, 32 

-sm - -sm + - 



8r./..N 8,^,,,, 40^,,,, 8, 224^,,,, 176 



(7.23) 



R'-^F' 



2 _ 5 
7 ~ 6 



(7.24) 
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8 Appendix B: Finite and Infinite Sums 



In the following, we list several classes of sums, which were used to derive the results in the 
present paper beyond well-known results for harmonic sums. Other relations can be found 
in [27,28,34-36] and were used in the present calculation. Here N, L, A denote arbitrary integers, 
a is a complex number, and B{a, b) is Euler's Beta- function. 

8.1 Sums involving Beta- Functions 



B{N,i) 



i=l 

oo 

1=1 



B{N,t) 



i=l 
oo 



i + 2 



^B{N,i) 



i=l 



E 

i=l 
oo 

E 

i=l 
oo 

E 

i=l 



i + 3 



B{N,t) 
B{N,i) 



B{N,i) 
N + i 



^ 1 + N + i 

2 + N + i 

B{N,t) 

3 + N + i 

BiN,t) 

4 + N + i 



1=1 

oo 

E 

1=1 

oo 

E 

i=l 
oo 



C2-S2{N-1) , 



1 + 



S2{N-1)-C2 {N-1), 



5-2N 



+ 



6N^ - 33N + 49 
36 

+ k(iv-i)-C2 



{N - 1){N - 2){N - 3) 



= C3-C2Si{N-l) + Si,2{N-l) , 

= la + C2^i,i(A^ - 1) - ^i,i,2(A^ - 1) - C3^i(A^ - 1) , 
o 

_ 1 



i=l 



7V4 + 4iv3 + 7iv2 + 6iV + 4 
m{N+iy{N + 2f ' 

+ 9N^ + SAN^ + 69N^ + S5N'^ + 66N + 36 
A^2(7V + l)2(7V + 2)2(Ar + 3)2 

m) 

m{N + l)2(iV + 2)2(iV + 3)2(iV + 4)2 ' 
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P{N) = Ar^ + 16A^^ + 110A^^ + 424Ar5 + l013A^^ + 1576A^=^ 



E 

i=l 
oo 

E 

i=l 
oo 

E 



B{N + l.i) 
N + i 

B{N + 2,i) 
N + i 

B{N+l,i) 



E 



E 



E 



^ (2 + Ar + i)2 



+1660A^^ + 1200A^ + 576 , 



1)^ 2S_2{N) + C2 



1 ' 



(-1)^ [Ca + 5i(A^)C2 + 25i,_2(A^) + 5-3(A^)] 

(-1)^ 



N{N + 1) 



25_2(iV)+C2 



2(_i)^25'_,(A- + 2) + C2 



+ 



+ 



N-1 
N{N + iy ' 

A-^ + A- + 1 



A^(A^ + l)(Ar + 2) N{N + lf{N + 2)^ ' 



■1) 



AT 



C3 + ^i(Ar + 1)C2 - C2 + 2S,,_2{N + 1) 



A^(A^ + 1) 
-2S_2{N + 1) + S_3{N + 1) 



Sums involving Beta-Functions and Harmonic Sums 

00 

J2B{N,^)Si{i) = C2-^2(A^-2), 
B{N,i)Si(i) 



i=l 
00 

E 

i=l 



E 

i=l 

00 

E 

i=l 
00 

E 

i=l 

00 

E 

i=l 



B{N,i)Si{i) 



N + i 

B{N,t)Si{t) 
N + i + 1 

N + i + 2 



2C3 + Si(N - 1)S2(N - 1) - C2Si(N - 1) 
-S2,i{N-l) , 

^Cl + C2^i,i(^-l)-2C3^i(^-l) 
-S,,,,2{N - 1) + S,,s{N - 1) , 



B{N,t)S,{t) C2-S2{N-1) 



N 

C2-S2{N-1) 



+ 



1 



N + 1 ' m{N + i) ' 



N + 2 



N^{N + lf 
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E 

i=l 



B{N + l,i)Si{z) 
N + i 



i=l 

1=1 



i=l 



B(A'.;)Si(A' + , - 1) 



E 

i=l 
oo 

E 

i=l 



N + i 

B{N,i)Si{N + i) 
N + l + i 



E 

i=l 



E 

i=l 



B{N,z)Si{N + z) 
N + 2 + i 



^ {N + i + iy 



i=l 



E 

i=l 



D{N + l.i)S,{N + i) 
N + i 



N 



N 



(2 



+2S^MN) - ^ + C2Si{N) 
Si{N- 1) ^ 2A^2 _ 2iv+ 1 



A^- 1 



m{N-iy ' 



2C3-253(iV) + 5i(iV) [C2-^2(iV) 

1 



^i(iV) 



+ 



2C3 - 253(iV - 1) + S,{N - 1) C2 - ^2(iV - 1) 



B{N,i)Si{N + i) Si{N-l) 2 



N^ 

N^ + N + 1 

m{N + if 
1 



+ Ar3 ' 



Si{N) - 



■1) 



AT 



N{N + 1) 



2S.2{N)+C2 



B(N,i)Si(N + i- 1) _ N'^ + N+1 
N+l + i ~ N'^{N + iy 

2N + 1 



Si{N) - 



(-1) 



AT 



A^(A^ + 1) 



2S_2{N)+C2 



N%N + iy ' 

N^ + 2N'^ + 5N + 2 
Ar3(Ar + 1)2(2 + AT) 

+ 4Ar3 + 77V2 + 6A^ + 4 



+ 



Ar2(7V + l)2(7V + 2)2 
2(-l)^ 



^i(iV) 



A^(A^ + l)(A^ + 2) 

(-ir 



2S_2{N) + C2 



25_2,i(7V+l)-25i,_2(A^ + l) 



Ar(A^ + l) 

+45_2(Ar + 1) + 2C2 + Cs - (2-51 (AT + 1) 

Si{N+l) 1 
^Ar(Ar+l)2 + N{N + 1)^ ' 

(-l)^[25_2,i(iV) + 5_3(7V)+2C3 



25 



E 

i=l 



B{N,i)Si{N + i) 



42 



i=l 



B{N,i)Si{N + i- 1) 



+ C2 



SijN) 
N 



- 2^1,1 (iV) 



+S2,2{N) + 2Si,s{N) + Si{N)S,,2{N) - 2 
S,{N)S2{N) 



SsjN) 
N 



N 



(8.34) 



la - CsSliN - 1) - 2C2-Si,i(iV - 1) + S2,2{N - 1) 
+25i,3(iV -l) + Si(N- l)Si,2{N - 1) , (8.35) 



i=l 
00 

E 

i=l 



B(N,t)Sf{i) 



E 

1=1 



E 

i=l 

oc 



BiN,i)S2ii) 



B{N,t)S2{i) 
N + i 



^ 1 + N + i 
1=1 

00 

J2B{N,t)S,,,{N + t) 



i=l 



1=1 



E 

i=l 



B{N,i)Si^i{N + i) 



3C3 - C2Si{N - 2) + 5i,2(iV - 2) - 2S'3(iV - 2) , (8.36) 



i=l 

00 

^S(7V,z)52(i) = -5i(iV-2)C2 + C3 + ^i,2(A^-2) , 



+ ^C2' - 3C3^i(A^ - 1) + 1(2 [sun - 1) 

-S2{N - 1)] - ,Si,i,2(A^ - 1) + -52,2(A^ - 1) 

+2S,(N - l)Ss(N - 1) - 253,i(7V - 1) , 
S^^,{N - 1) - 5i,i,2(A^ - 1) + C2^i,i(A^ - 1) 

-C2^2(A^-l)-C3^l(A^-l) + ^C2', 

-S,{N -1X2 + (3 + Si,2{N-l) 
N 

-S,{N)C2 + C3 + ^i,2(A^) ^2(iV) - C2 



1 



J2B{N,^)S^,^{N + ^-l) = 



N+1 

gi,i(iV-l) S^{N-l) 

N-1 (AT - 1)2 + {N-iy 

, ^i(iV-l) , 2 

ATS ' 

^i,i(Ar-l) ^ 5i(Ar-l) ^ 1 



A^- 1 



(A^-l)2 (A^-l)3 

Si,i{N) C2 - 52(A^)] + Id - ^S4N) 



+2S,{N)[C3-S3{N)\ + 



Si{N) ^Si,,{N) , 1 



(8.37) 



(8.38) 

(8.39) 
(8.40) 
(8.41) 



(8.42) 
(8.43) 



7V3 
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i=l 



B{N,t)Si,i{N + i-l) 



Si,i{N - 1) C2 - ^2(A^ - 1) + - ^S,{N - 1) 



+2S,{N-1) (^-S3{N-1) 



i=l 

00 

1=1 



B{N,i)S,,,iN + i) 
N + i 

B{N,i)Si^,{N + i) 
N + i + 1 



1 

TV 



S,,,{N) ^ S,{N) ^ 1 



N 



S^{N) S,,^{N){N^ + N+l) 1 
ATS + Ar2(Ar + l)2 



N 



N{N + 1) 



2C3 + 5_3(iV)+25_2,l(iV) 



i=l 



^(7V,i),gi,i(7V + i-l) 



,5i(A^ ~ ^) + - 1)(A^2 + + 1) 



A^2(7V + l) 
(-1)^ 



C2 + 2C3 + 25_2(7V-1) 



A^(A^ + 1) 
+S-3{N -l) + 2S.2,i{N -I) 



00 

^B(Af,i)Si(i)Si(Af + i) = S,(A'-l)[c2-S2(Af-l)] +2[C3-53(A'-1)] 



E 

1=1 



S(Ar,i)5i(i)5i(Ar + i) 



C2 _ ^2(A^-l) >gi(iV-l) 2 
AT AT ^ (AT - 1)2 + (iV- 1)3 



17 

To 



,C3 

'N 



(8.45) 



(8.46) 



(8.47) 



(8.48) 



(8.49) 



1 ^£i(^^2S„(iV-i) 



LA^2 



N 



^ S2{N - 1) _ 53(iV - 1) ^ _^i,2(A^ - 1) 



Ar2 



N 



N 



+S,{N - l)(^Ss{N - 1) + S,,2{N - 1)) 

+>52(Ar - 1)2 - 52,2(A^ - 1) - 2>53,i(A^ - 1) , (8.50) 



E 

1=1 



E(Ar,^)^i(^)^i(Ar + i) 
AT + i 



E 



B{N.i}S,{i}S,{N + {) 
1 + N + i 



j^{s,{N) [C2 - S2{N)'\ + 2[C3 - 53(A^)] 



^Si(AO 2 ^ 
7V2 ^ iV3 J ' 

'_1)^+1 



Ar(Ar + l) 

, C2^l(A^) 



Cs + C2^i(A^) + ^-3(^) + 25i,_2(Ar) 



A^ + 1 

2Ss{N) 

N + 1 ' m ' 



+ S,{N) 
2 

+ 



_]__S2{N)_ 
N + 1 



+ 



2C3 



A^+ 1 



(8.51) 



(8.52) 
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T{N + k) B{k + i,N) Si{i) 17 



^^r(7V)r(/. + i) k 



10 



Cl-2CsSi{N-l) + S2,2{N-l) 



+Si{N - l)S2,i{N - 1) - 2S2,iAN - 1)(8.53) 



8.3 Sums involving Beta- Functions and Harmonic Sums with two 
Free Indexes 



^B{N -a,i + a) = B{1 + a, N - 1 - a) , 



(8.54) 



1=1 



E 



B{N, i + k) 



d 



dN 



B{k,N) 



= B{k,N) Si{k + N-l)-Si{N-l) , (8.55) 



i=l 



B{k, N) 2Si^i{k + N-1)- S2{N - 1) 



-S^{N + k-l)S,{N-l)\ , 



(8.56) 



y ^^^''^^^ ^i(i + A; - 1) = S(A;, AT) Si{k + N - l)Si{k - 1) - 52(A^ - 1) 



-5i(A;-l)5i(iV-l) + C2 



(8.57) 



8.4 Sums involving Binomials 



A;=0 



-1)' 



B{N,A) 



{k + Ay 



[s,{N + A-l)-S^{A-l)Y 



-S2{A-l) + S2iN + A-l) 



- iim{m.)-i}=-5,(iv-i), 

J.— 1 V / 



k=l 
N-l 

E 

k=l 
N-l 



A; r ' k^ 



k=l 



k r k^ 



E(^:^)(-l)^ 



N-l 



k=l 



1 



{k + iy 



= -5i,i(iV-l), 



-5i,i,i(iV-l) , 



N 



(8.58) 
(8.59) 
(8.60) 
(8.61) 
(8.62) 
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A; ^ (A; + 2)3 Ar(Ar + 1) (AT + 1)2 (AT + 1)3 8' ^' ^ 



N-l 



J^(^-lV_l)fc 1 - 25i,i(iV) (5 + 3iV)^i(iV) 



ik=l 



(A; + 3)3 iV(Ar + i)(iV + 2) (A^ + l)2(Ar + 2)2 

^(Ar + l)3(Ar + 2)3 ~ 27 ' ^ ^ 

^[N^l\ , 1 ^ >gi,i,i(iV) S^,,{N) S,{N) 1 

J (A; + 2)4 Ar(iV + l) (iV + l)2 (A^+l)3 {N + l^ 

1 

~16 ' 

"^"^ ^AT-lX 1 _ 25i,i,i(Ar) (5 + 3Ar)5i,i(Ar) 



ik=l 



E -r (-1) 



A; ' (A; + 3)4 Ar(Ar + l)(Ar + 2) (AT + l)2(Ar + 2)2 

(Ar3- SAT -9)^1 (AT) 



+ 



(A^ + l)3(A^ + 2)3 



8.5 Sums involving Binomials and Harmonic Sums 

We set Sa{0) ■= 0. Therefore the sums below may begin at A; = 0. 

N-l 



k=l 

N-l 



k=0 
N-l 



k=0 

N-l 



N-l\. ^,fc5i(A; + 2) S2{N) 



k=0 



(8.66) 



-17-21A^-2Ar2 + 6A^3^2iV4 1 
^ (Ar + l)4(A^ + 2)4 81 ■ ^ ' ' 



X^rr)^-^)'^ - lim{5(iV,e)[5,(Ar + .-l)-5,(Ar-l)]} 

= C2 - S2(N - 1) , (8.68) 

E^r^V-l)'^ - -S.AN-1), (8.69, 



^\ k )^ ' {k + 2f N{N + l) (Ar + 1)3' ^ 
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k=0 



N-l 



E 

fe=i 

E 

k=l 
N-l 

E 

fe=l 

N-l 

E 

k=l 



N-l 
k 

N-l 
k 

N-l 
k 

N-l 
k 



k + 2 



-1) 
-1) 
-1) 
-1) 



kS2{k) 



k 

k S2{k) 

k+1 

2 + k 

3 + k 



N-l 



k=l 

N-l 



^.N-l\ , S,ik) 



k=l 



{2 + ky 



yV^^-iy ..k S^jk) 



N-l 



fe=0 



k 



N-l 

E 

k=0 

N-l 



k + A 



N - 1 
k 



■1)' 



■ Si,i{k) 
k 



k=i ^ ^ 



N-l 

E 

k=l 



N-l 
k 



k+1 



;bi[N) + ■ 



N^(N + iy 



= -^2,l(iV-l) , 



7V2 



Ar3(7V + l)3 



AT 



1- AT 



2^1,1 (iV) 



SiiN) 



1 



N+1 ' 
+ 4Ar - 4 



(8.73) 



(8.74) 



(8.75) 



(8.76) 



A^(A^ + l)(A^ + 2) 2{N + 2)m 
N + 11 



Si{N) 



4:{N +1){N + 2) ' 
S2,i{N) -Si,i,i(7V) 



N 



N 



S2,i{N) S,,i,,{N) , Si,,{N) N + 2 



+ 



N{N + 1) N{N + 1) {N + iy N{N + 1) 
2N + S 



2^i,i,i(^) 



^(7V + 1)2 ' 

2^2,i(iV) 

N{N + 1){N + 2) N{N + 1){N + 2) 

{3N + 5)Si,i{N) _A0 + 38N + 9N^ + N^ 
^(A^+ l)2(iV + 2)2 ~ AN{N + 1){N + 2y 
59 + 66iV + ISiV^ + iV3 



(8.77) 
(8.78) 

Si{N) 
(8.79) 



5i(A^) 



4(iV + l)2(Ar + 2)2 
B{N, A) 



(8.80) 



|5i(A^ + A-l)-,5i(Ar-l)y 
+,S2(A^-1 + A)-52(A^-1)] , 

Cs - ^3(iV - 1) , 

-S2,i,i(N-l) , 
Si,iAN-l) 



N 
30 



(8.81) 

(8.82) 
(8.83) 
(8.84) 



N-l 

E 

ik=0 



N-l 
k=l 



N-l 
k 



N-l 
k 



' 3 + k 



Er.:^i(-:) 



fe=0 

N-l 



, .fc >gi,i,i(A: + A) 
^ ' k^A 

A\^{k) 



k 



fc=0 

Af-l 



k 

.^2,1 (A:) 



A; 



5iii(iV) l-AT^ ^^^^ 5i(Ar) 



A^(A^ + 1) 
1 



Ar+ 1 ' 

25i,i,i(A^) _ iV^ + 4iV - 4 
'iV(iV + l)(Ar + 2) ~ 2(Ar + 2)iV2 

^i(Ar)(Ar + 10) iV+19 



(8.85) 



47V(Ar + 2) 8(7V + l)(Ar + 2) ' 

1 d3 



5(iV,^) , 



7 
10 



a - S2AN - 1) , 



fe=0 ^ ^ k 



A;2 



(8.86) 
(8.87) 
(8.88) 
(8.89) 



= -^C2 + C3^i(A^ - 1) - 2C25i,i(iV - 1) + 5i,3(A^ - 1) + 2S^,,,2{N - 1) . (8.90) 
8.6 Sums involving Binomials with Two Free Indexes 



L+l 

E 



fc=0 



L+l 



k jN-L+k 



L+l 



E 

fe=0 
L+l 

E 

fe=0 



L + l 
A; 

L + l 
k 



■1) 



^ 5^(7V-L + fe) 
N-L + k 



, S,{N-L + k) 
1+N-L+k 



B{N-L,L + 2) , 

B(N-L,L + 2) ^Si{N + L) 
-S^(N-L-1)\ , 



B{N - L,L + 2) Si{N +1) - Si{L + l) , (8.93) 



(8.91) 



(8.92) 



B{N -L + l,L + 2) Si{N-L)- Si{L + 1) 
d 



dL 



B{N - L + l,L + 2) , 



(8.94) 



31 



E 

k=l 



B{k + e/2,N + l) 
N + k 



[-1 



25_2(7V) + C2 



-2 



-Ca + Q2Si{N) + 25i,_2(iV) - 25_2,i(7V) 
2 



+ - C3^i(iV) + C2S,,,{N) + 2{5i,i,_2(iV) 

+^-2,l,l(A^)-^l,-2,l(A^)}]+0(£3) , 



.95) 
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Figure 1: Feynman rules for the operator insertion ® to O(a^), cf. [21]. A denotes a light-like 
vector, A. A = 0. 




Figure 2: The Feynman diagrams contribution to the operator matrix element Aqg at 0{as). Weavy 
lines denote gluons, and the full arrow lines are the heavy quark lines. The Feynman rules for the 
operator insertions are given in Figure 1. 
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Figure 3: The diagrams contributing to the operator matrix element Aqg at O(a^). Weavy lines 
denote gluons, dashed lines ghosts, and the full arrow lines are the heavy quark lines. The Feynman 
rules for the operator insertions are given in Figure 1. 



37 



Figure 4: The diagrams contributing to the operator matrix element at O(a^). Weavy lines 
denote gluons, the thick full arrow lines are the heavy quark lines, and the thin full lines are light 
quark lines. The Feynman rules for the operator insertions are given in Figure 1. 




Figure 5: The diagrams contributing to the operator matrix element Aq^^ at O(a^). Weavy lines 
denote gluons, the full arrow lines are the heavy quark lines, and the thin full lines are light quark 
lines. The Feynman rules for the operator insertions are given in Figure 1. 
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